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This paper introduces a significant special situation where the noise is a 
collection of D-plane interference signals and the correlated noise of D+1 is 
less than the number of array components. An optimal beamforming 
processor based on the minimum variance distortionless response (MVDR) 
generates and combines appropriate statistics for the D+1 model. Instead of 
the original space of the N-dimensional problem, the interference signal 
subspace is reduced to D+1. Typical antenna arrays in many modem 
communication networks absorb waves generated from multiple point 
sources. An analytical formula was derived to improve the signal to 
interference and noise ratio (SINR) obtained from the steering errors of the 
two beamformers. The proposed MVDR processor-based beamforming does 
not enforce general constraints. Therefore, it can also be used in systems 
where the steering vector is compromised by gain. Simulation results show 
that the output of the proposed beamformer based on the MVDR processor is 
usually close to the ideal state within a wide range of signal-to-noise ratio 
and signal-to-interference ratio. The MVDR processor-based beamformer 
has been experimentally evaluated. The proposed processor-based MVDR 
system significantly improves performance for large interference white noise 
ratio (INR) in the sidelobe region and provide an appropriate beam pattern. 
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1. INTRODUCTION 


Beamformers are commonly used in a variety of applications of signal processing, such as voice 
amplification, radar, and wireless communication. In particular, one of the fundamental problems in array 
processing is the problem of measuring the signals of several transmitters using an antenna array. The Capon 
beamformer [1] has gained popularity in this context due to its high ability to refuse interference. 
Nevertheless, Capon beamformers also known as Van Trees minimum power distortionless response 
(MPDR) beamformers are still prone to steering errors [2]. Ward, Cox, and Kogon based on diagonal loading 
effect provides a comparison of several robust MPDR beamformers to correct model errors [3]. Wax et al. 
[4] performed a theoretical score reading of the signal-to-interference-plus-noise ratio (SINR), which is 
generalized (that is, not necessarily a diagonal) covariance matrix and the existence of random steering vector 
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errors. Interestingly, diagonal loads have proven to be the worst-case solution proposed recently 
[5]-[7]. The beamformer is built in the latter references to minimize the output power subject to the 
restriction that the response of the beamformer is above some level for all the steering vectors in an ellipsoid 
centered around the nominal or presumed value steering vector. In the presence of partially and fully 
correlated interfering sources, Reddy et al. [8] first study the signal cancelation and interference rejection 
effects of the optimum (constrained minimum squares or minimum variance) beamformer. 

Gu and Leshem, proposed an effective adaptive beamforming algorithm [9]. This algorithm is 
robust not only to the covariance matrix uncertainty but also to steering vector mismatch. Using direction of 
arrival (DOA) knowledge to estimate the covariance of interference and noise that can reconstruct a matrix 
from a spatial spectral distribution. In the presence of unknown signals, a new robust and adaptive 
beamforming method was developed by Vorobyov et al. [10] to steer vector mismatch [11]. This method is 
based on worst-case performance optimization. The natural formulation of this adaptive beamforming 
problem involves minimizing a quadratic function with an infinite number of non-convex quadratic 
constraints. The main classification is for a system dedicated to dealing with covariance lattice examples, 
because apart from the clamor covariance framework, precise obstacles are usually not accessible in 
near-Earth applications. The best known of this class is called the diagonal loading method [12], [13] where a 
scaled identity matrix is added to the sample covariance matrix. However, choosing the right diagonal load 
factor for different situations can be a daunting task. Recently, the shrinkage estimation method [14], [15] 
automatically measures the diagonal load factor without defining user parameters. Unfortunately, this only 
provides an approximation of the theoretical received signal covariance matrix, not the required interference 
plus noise covariance matrix. The second method classifies the signal’s expected steering vector explicitly, 
because it is not easy to obtain accurate steering vector information. In this case, the adaptive beamforming 
technique based on worst-case performance optimization [16], [17] explicitly uses an uncertainty array of 
signal control vectors. In fact, the mismatch vectors and their boundary norms are unknown. In this case, the 
most recent method is to iteratively estimate the actual steering vector [18], [19]. 

If the array variants are fully known, there is no difference between using the acquired signal 
covariance matrix MPDR or sample matrix inversion (SMI) and using minimum variance distortionless 
response (MVDR) [20]. However, if sequence variants are not completely clear, we recommend using the 
MVDR method [21], [22]. On the other hand, most recently developed robust technologies use a 
generalization of MPDR technology. Khabbazibasmenj et al. [23], [24] and Mallipeddi et al. [25] recent 
efforts have been made to isolate the effects of interference. This form of matrix replaces the interference and 
noise covariance matrices in two published articles. Reference [26] changed the optimal signal control vector 
to only be used as part of the semi-optimal optimization. This increases the signal loss of the noise output. In 
[27], the integral is replaced by a discrete DOA, which is determined by minimizing the Capon spectrum 
estimator. The fact that interference power estimation is not included in the construction of the covariance 
matrix undermines methods used for high dynamic range applications, such as radio astronomy imaging of 
diffuse sources. The response is based on the reconstruction of the interference plus noise covariance matrix 
and the steering vector approximation to evaluate the optimal beamformer for discrete interference [28]. As 
mentioned above, in order to arrange multiple errors, the MVDR beamformer is more flexible than the 
MPDR beamformer. New design ideas for multiple-input-multiple-output (MIMO) beamforming antenna 
arrays for compact and thin handheld devices reported by Liu et al. [29]. A printed array clutch broadband 
beamforming with four 2.4 GHz multiple input multiple output chip antennas (CA) for 5G communication 
has been investigated by Hansen et al. [30]. 

Therefore, under the constraint that the modified steering vector does not converge to interference, 
the assumed steering vector of the signal is subsequently corrected to optimize the output power of the 
beamformer. By combining the interference-noising covariance matrix and the adaptive beamformer, a 
reconstruction based on steering vector estimation can be obtained. These also do not require general 
restrictions on steering vectors [31], [32]. 

Therefore, beamforming based on MVDR processors can be extended to more situations than 
before. Numerical examples show that for low signal-to-noise ratio (SNR) and high SNR performance, using 
the previously proposed robust beamformer, the performance of the proposed beamforming algorithm is 
almost always the best [33], [34]. Combining the reconstruction of the interference plus noise covariance 
matrix with the enhanced estimation of the required signal steering vector can overcome the problem of 
self-cancellation of the required signal at high SNR, while maintaining good performance at low SNR. This 
improves performance compared to existing methods and has a wide range of signal-to-interference and noise 
ratios [35]. In operation, multiple signals from a point source are received by the array. Interference sources 
and non-coherent sources form a noise correlation matrix. The results show that the MVDR beamformer is 
still more robust in the presence of uncertain interfering signal steering vectors [36]—[38]. 
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2. SIGNAL COPY AND OPTIMUM DETECTION PROBLEM FORMULATION 

If, as is typically the case, the noise and the desired underlying signal are uncorrelated, then the 
variance of the reported signal is the sum of the desired signal and noise variances. The MVDR solution is 
therefore aimed at reducing this number, thereby mitigating the noise effect. A solitary plane wave signal 
associated to the array manifold v(w: K,)and a signal that interferes with the array manifold v(w: K1), plus 
the white noise component of individually sensor at the aĝ spectral height. Our main goal to find the best 
undistorted filter, the best typical beamformer array gain, and the best beamformer beam pattern. The block 
diagram of MDVR processor is shown in Figure 1. A case of plane wave interference reveals several 
outcomes that can be generalized to many interfering signals. 


S,(@) = 021 + M,(w)v(@: K,)v" (w: K1) (1) 


In (1) M; (w) is the spectrum of interfering signal S, = 021 + M,v,v!! can be obtained by suppressing w and 
K. Applying the matrix inversion lemma. 


— i = Mı H 
Sn e on, | o2, +NMı Vivy l (2) 


Expected noise at each sensors estimated as follows: 


S, = (02 + My) (3) 


X(w) A(w)v" (w: K,)S,+(w) Vw 


Figure 1. Block diagram of MVDR processor 


We ay > WE (w) =A (w)v" (w: K,)Sz1(w) 
(4) 


The ratio of the signal to the noise spectrum at the output of the undistorted filter relative to the 
input ratio gives the array gain at a specific frequency w. Replace the value of (2). In (4), obtain Wyar as 
shown in (5). 


H —_,,H ^H Mı H 
Wnmvar =Wo = z2 Us |i B ViVi | (5) 
w 


vg V1 

N 
represents predictable beam pattern expected as Ks which is the signal wavenumber and estimated at Kis the 
interfering wave number. Then (5) can be reduce as follows: 


A 
Defining spatial correlation coefficient ps, > among the anticipated signal to interference ratio. ps, is 


H AN p NM, vt) 
Wo oZ LN Psi o2+NM,° N (6) 
Where, A is the normalizing coefficient given as: 
^= | 1 yH {1 M y vi} v E = |= Nfi NM eye (7) 
e 03, > o +NM1 ae 2 a 02, oZ +NM1 Psı 


Optimum MVDR processor of single plane-wave interferer for general and high interference case is 
shown in Figure 2(a). We see that the beamformer generates two sufficient spatial statistics, one that 
corresponds to the anticipated signal and one that corresponds interfering signal, and then optimally 
combines them. Considering that each of the space processors is, respectively, a traditional beamformer 
pointing to the signal and interferer. Then we deduct a fraction of the estimate that corresponds to the two 
signals spatial correlation. In case NM, > aå then WË can be derived as [39]: 
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— 1 
Wé = vd [I -vivn tof] = 4 vl Pe (8) 


where p1 is the interference projection matrix on the orthogonal subspace. Figure 2(b) displays the 


corresponding device. The beamformer puts on the interfering signal a complete zero. Substituting 
1 NM ees z 

Sn = (o2 + M,) and [an {1 - a losi? | into Ao(w:ks) = Sn (w)v” (w: ks)S7 (wv (w: ks) to 

obtain optimum array gain. 


Ao =n(1+) (1 - ilps?) (9) 


Xto) => Y(v) 


(b) 


Figure 2. Optimum MVDR processors: Single plane-wave interferer general and high interference case (a) 
optimum MVDR processor; single plane-wave interferer general case in correlated signal and interferences 
and (b) optimum MVDR processor; single plane-wave high interferer 


2 
Ao = NG + of) [rÉ (10) 


2 
1+Nof 


It is found that the optimal gain of the array depends on øf, INR and sy, that is, the spatial 
correlation between the signal and interference. Computer simulated first 101 optimum array gain for 
10-element array are as follows: 1.009, 11.099, 21.189, 31.279, 41.369, 51.459, 61.548, 71.638, 81.728, 
91.818, 101.908, 111.998, 122.088, 132.178, 142.268, 152.358, 162.448, 172.537, 182.627, 192.717, 
202.807, 212.897, 222.987, 233.077, 243.167, 253.257, 263.347, 273.437, 283.526, 293.616, 303.706, 
313.796, 323.886, 333.976, 344.066, 354.156, 364.246, 374.336, 384.426, 394.515, 404.605, 414.695, 
424.785, 434.875, 444.965, 455.055, 465.145, 475.235, 485.325, 495.415, 505.504, 515.594, 525.684, 
535.774, 545.864, 555.954, 566.044, 576.134, 586.224, 596.314, 606.404, 616.494, 626.583636.673 646.763 
656.853 666.943 677.033 687.123 697.213 707.303 717.393 727.483 737.572 747.662 757.752 767.842 
777.932 788.022 798.112 808.202 818.292 828.382 838.472 848.561 858.651 868.741 878.831 888.921 
899.011 909.101 919.191 929.281 939.371 949.461 959.550 969.640 979.730 989.820 999.910 1010.000. 
Figure 3 represents array gain versus {1 — |ps1|7}: 10-element array for INR = 0 dB, INR = 20 dB and 
INR = —20 dB in correlated signal and interferences. The improvement of the array gain has been observed 
24dB for INR = 20 dB and 5.343dB for INR = 0 dB at |ps,|? = 0.5 half power beam width (HPBW) line. 
The improvement of array gain is marginally high for INR = 20 dB at |ps5,|* = 0.5 and |p.,|* = 0.95 line 
compared to INR = OdB. Using the standard beam pattern ps, = B,(k,,k,;), array geometry must be 
incorporated in ps1 so that these tests are applicable to an arbitrary array. HPBW lines at |ps,|? = 0.5 and 
|Ps1|* = 0.95 is shown in Figure 3. The interferer is inside the main lobes HPBW region for |ps1|* = 0.5 to 
the left of this axis. Substituting Sp = 021 +M,v,v! into A.(w:k,) to get the conventional array gain 


A, (w: ks) = N(1+07) 


1407 Nlpsil?" 
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35 T 
—— Optimum for INR -20dB 


= = Conventional for INR -20dB 
= = Conventional for INR 0dB 
= Optimum for INR 0 dB 

= = Conventional for INR 20dB 
= Optimum for INR 20dB 
—— HPBW line at |p, 1 |?=0.5 


30 


N 
a 


= = HPBW line at |p, 1 |7=0.95 


N 
o 


~- — 


= 
a 
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INR=0 dB 


Array gain (dB) 
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= 
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= 
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- 
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Figure 3. Array gain versus {1 — |ps1|7}: 10-element array: For= 0 dB, INR = 20 dB and INR = —20 dB 
in correlated signal and interferences 


3. PATTERN PROBLEM ANALYSIS 

Noticed that through the standard beam pattern in the array geometry is encoded in ps, so that these 
effects are applicable to an arbitrary array. The interferer is inside the main lobes half power beam width 
(HPBW) region for |ps1|? to the left of this line. When we look at beam patterns in the examples, we find 
that this can cause pattern problems. Considering the number of elements N is very large then 
No?Z[1 — |ps1|*] > 1 and optimum array gain Ap becomes as in (11). 


Ap = N(1 + 07)(1 — |ps1 I”) (11) 


Computer simulated first 100 optimum array gain for traditional system for 10-element array are as 
follows: 0.010, 0.020, 0.030, 0.040, 0.051, 0.061, 0.071, 0.081, 0.091, 0.101, 0.111, 0.121, 0.131, 0.141, 
0.152, 0.162, 0.172, 0.182, 0.192, 0.202, 0.212, 0.222, 0.232, 0.242, 0.253, 0.263, 0.273, 0.283, 0.293, 0.303, 
0.313, 0.323, 0.333, 0.343, 0.354, 0.364, 0.374, 0.384, 0.394, 0.404, 0.414, 0.424, 0.434, 0.444, 0.455, 0.465, 
0.475, 0.485, 0.495, 0.505, 0.515, 0.525, 0.535, 0.545, 0.556, 0.566, 0.576, 0.586, 0.596, 0.606, 0.616, 0.626, 
0.6360.646 0.657 0.667 0.677 0.687 0.697 0.707 0.717 0.727 0.737 0.747 0.758 0.768 0.778 0.788 0.798 
0.808 0.818 0.828 0.838 0.848 0.859 0.869 0.879 0.889 0.899 0.909 0.919 0.929 0.939 0.949 0.960 0.970 
0.980 0.990 1.000 1.010. Therefore, plotting Ag/N as shown in Figure 4(a). The beam pattern can be 


: ae AN [vë NM, vf]. es r : 
obtained by substituting Wo’ = -z É = Psi zonm a] in the probability density of X (wm), given F (wm) 
w w 1 


is for the optimum array represented in mathematical form as follows: 


PX (Wm) |F (wm) C) 
= cyexp[—{X" (wm) — F* (wm) V” (®m)}Sn* (@m)[X (wm) 


— F (wmv (wm)]] (12) 
The result is 
N No? ps 
Bolo: ks) = S [psa — SEP pra (13) 


where A is the normalizing coefficient and suppressing the w in right side of (13), 
1 
Psa = —v"(w:ks)v(w:k) = B,(k: ks) (14) 


which is the traditional beam pattern for a k,-pointed delayed and summed of beamformer and that can be 
expressed mathematically as follows: 
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Pra = 5, V” (w: ky) v(w: k) = Belk: ki) (15) 


The (15) is only the typical beam pattern for a k;-spot beamformer. Its optimum array pattern can be further 
interpreted as the traditional beam pattern pointing to the desired signal minus a constant time pointing to the 
interference by the beam pattern. 


A(@).N 


B(w, k: ks) = 


2 
B.(w, k: ks) — aa B, (kı: ks). Be(k: ky) (16) 


o ' 


The value of the multiplication factor depends on two elements. The first is the k, value of the 
standard beam pattern pointing to k,. If the interference is zero in legacy mode, it has no effect and does not 
need to be subtracted. If the source of interference is in the traditional sidelobe, the corresponding sum must 
be deleted. The second term depends on a value No? near unity. The effect of combining these two items is 
to create partial zeros in the kų beam pattern. The beam pattern value of kı is mathematically expressed as: 


Alw) 


ot” 


Bo (key: ks) = “DNB. Ck: ko) (25) (17) 


1+No} 


As No? > © one can get the optimum array pattern null in the direction of interferer [40]. 


4. ANALYSIS OF THE MVDR PROCESSOR BASED BEAMFORMER 
Consider the generic linear array of 10-elements along the z-axis. From Ay = N(1 + øf) 
(1 — |ps,|?) the standard beam pattern which is designed for an array that is pointed out at broadside u, = 0 
and u = cos@ is expressed as B,(u). The region and their range of lobes is shown in Table 1. 
. (Nru 
Ba) = 22 


N sin(—* 


4 


(18) 


YS 


Table 1. Region and their range of lobe 


Region Range 
Sidelobe Region 0.2 < |u| < 1.0 
Outer Main Lobe 0.045 < |u| < 0.2 

HPBW Region 0 < |u| < 0.045 


Let us consider that interference is arriving from u; which is shown in Figure 4(b) illustrates the plot 
of Ap as a function of u;/BWyy, Ao/A, and u; observed its extreme at the HPBW. It has been observed that 
the interferer is occurring outside the main lobe that is u; > 0.5 BWyy and the optimum array gain is at 
N.INR GB. The traditional gain in the array has the similar performance as the u-space. 


——A, at INR 0dB 
= = A, IA, at0dB 
— A, at INR 20dB 
= = A/A, at 20dB 


— INR=20 dB 
25 = = INR=0 dB 
~~ INR=-20 dB 


see [pet |?=0.25 
= = [pt [=0.5 


(dB) 
dB 


o 


A_IN 


Figure 4. Correlated signal and interferences normalized optimum array gain (a) Normalized optimum array 
gain versus {1 — |p.,|7} and (b) A)/A,versus u;/BWyy ; IN = OdB and 20 dBin 
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The best antenna beam patterns of several values of uy; in the side lobe region of 10 dB IINR (of) 
are plotted in Figures 5(a) and 5(b). Figures 5(c) and 5(d) demonstration the best antenna beam pattern when 
the same u; value INR(a7) is 20dB in the sidelobe region. It can be seen that the antenna beam pattern is the 
same as the conventional antenna beam pattern except for the area close to zero. 

Figure 6 shows the optimal antenna beam pattern for the main lobe of external interference. 
Figure 6(a) of = 10 dB, u; = 0.18; Figure 6(b) of = 10 dB, u, = 0.09; Figure 6(c) of = 20 dB,u; = 0.18; 
Figure 6(d) øf = 20 dB, u; = 0.09. However, in the case of u,=0.09, these two effects are deceptive. (i) The 
main lobe moves from zero to the left and its height is greater than unity. (ii) The height of the right lobe is 
approximately 3 dB. Figure 7 illustrates the array gain for different gê values for N > 20 as a function of 
u,/BWyy.. Outside the main lobe region, only the square beam pattern is the array gain. Therefore, for 20 dB 
INR, performance may be reasonable. However, as the INR increases, the optimal array gain (Ay = N. INR) 
continues to increase, but the deterministic array gain remains the same. 


Beam pattern (dB) 

8 è o 
Beam pattern (dB) 
8 è ë o 


en 
ò 
in 
(= 
o 
èn 
a 
S 
in 
o 
o 
in 
- 


(a) (b) 


Beam pattern (dB) 
8 B è o 
Beam pattern (dB) 
ë & & o 


di 
b 
in 
o 
o 
in 
= 
a 
b 
in 
o 
o 
in 
- 


(c) (d) 


Figure 5. Optimum antenna beam pattern for side lobe interference (a) of = 10 dB, u; = 0.30, 
(b) of = 10 dB, u, = 0.50, (c) of = 20 dB, u; = 0.30 and (d) oê = 20 dB, u,; = 0.50 in correlated signal 
and interferences 


Beam pattern (dB) 
ù wm à P 
(= o (= (= o 
Beam pattern (dB) 
& wm + a 
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ò 
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o 
o 
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= 
o 
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& N = =- 
(= (= (= (= (= 

Beam pattern (dB) 
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o o (= iz (= 


1 05 0 05 1 1 05 05 1 
u u 
(c) (d) 


Figure 6. Optimum antenna beam pattern for outer main-lobe interference (a) of = 10 dB, u; = 0.18, 
(b) oÊ = 10 dB, u, = 0.09, (c) o7 = 20 dB,u,; = 0.18 and (d) of = 20 dB, u; = 0.09 in correlated signal 
and interferences 
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4.1. White noise plus spatially spread interference consideration 

Let us consider a 31-element standard linear array for MVDR beamformer. We assume that 
INR = 20 dB, INR = 10 dB, INR = 0 dB INR = —10 dB and INR = —20 dB are five interferers. There 
are two uncorrelated interfering signals about Y = 0 that are symmetrically located. Figure 7 illustrates the 
array gain for different øf values for N > 20 as a function of u;/BWyy. In the existence of white noise plus 
spatially spread interference, reflect the output of MVDR beamformer performance. Considering a regular 
linear array of 20-elements that is oriented vertically (along the z-axis). The snapshots of the frequency 
domain can be mathematically represented as: 


X(w) = vF (w) +N. (w) + Ww) (19) 


The interference comprises a white component with a spectral matrix 02] and a high-surface noise 
component with a spatial spectral matrix as: 


Sn,(@) = So(w) [sinc (ko Ap) + ja ap (sine (koAp) = cos(kgAp)}cos6,| (20) 


in this example, the surface noise is uncorrelated with the white noise. The sensor spacing half wavelength 
with respect to wc. The block diagram to obtain the optimum array gain is shown in Figure 8. 


25 T T 


INR=20 dB 


20 F 


Br INR=10 dB 


(dB) 


o 
a 
T 


A_IN 


INR=-10 dB 
INR=-20 dB 


0 0.5 1 1.5 2 2:5 
u, IBWyN 


Figure 7. Array gain for different øf values for N > 20 A/N (dB) as a function of u;/BWyyn 


Evaluating So (w: 0, Ø) Obtain the spectral matrix Compute S}! (w) = 
at the sensor location Sn, (w) [Sn (as) + J“ 
c 


N? 


Ao(@: ks) = v” (w: ks) pn'(w)v(w: ks) 
v4 (w;ks)Pnlw)v(w:ks) 


A,(w:k,) = 
elw: ks) By suppressing w and k, one can get 


Suppressing w 


i N . 2 - 
and k, to obtain A, = F PE for comparison Ao = V} Pr Ys 
sPn’s 


Figure 8. Block diagram to obtain the optimum array gain 


Computer generated 10x10 spectral matrix S,.(w) and S l@)= [Sn.(w) + VAR are shown in 
Table 2 and Table 3 respectively. Computer generated few values ofAo(w: ks) = v” (w: ks)pn (wv (w: ks) 
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N2 


by suppressing w and k,, Ap = vi pztv, and A,(w:k,) = LORIATA 


by suppressing w and k,, 


2 


N SOENE : 
Ac = —;—— for comparison is shown in Table 4. 
Vs Pnvs 


Table 2. Computer generated 10x10 spectral matrix S,.(w) 


Cl C2 C3 C4 CS C6 C7 C8 C9 C10 
R1 101.00 10.13 -2.53 1.13 -0.63 0.41 -0.28 0.21 -0.16 0.13 
R2 10.13 101.00 10.13 -2.53 1.13 -0.63 0.41 -0.28 0.21 -0.16 
R3 -2.53 10.13 101.00 10.13 -2.53 1.13 -0.63 0.41 -0.28 0.21 
R4 1.13 -2.53 10.13 101.00 10.13 -2.53 1.13 -0.63 0.41 -0.28 
R5 -0.63 1.13 -2.53 10.13 101.00 10.13 -2.53 1.13 -0.63 0.41 
R6 0.41 -0.63 1.13 -2.53 10.13 101.00 10.13 -2.53 1.13 -0.63 
R7 -0.28 0.41 -0.63 1.13 -2.53 10.13 101.00 10.13 -2.53 1.13 
R8 0.21 -0.28 0.41 -0.63 1.13 -2.53 10.13 101.00 10.13 -2.53 
R9 -0.16 0.21 -0.28 0.41 -0.63 1.13 -2.53 10.13 101.00 10.13 
R10 0.13 -0.16 0.21 -0.28 0.41 -0.63 1.13 -2.53 10.13 101.00 


Table 3. Computer generated 10x10 spectral matrix Sp + (w) = [Sn (w) + o2) 


Cl C2 C3 C4 C5 C6 C7 C8 C9 C10 
R1 0.0100 -0.0011 0.0004 -0.0002 0.0001 -0.0001 0.0001 0.0000 0.0000 0.0000 
R2 -0.0011 0.0101 -0.0011 0.0004 -0.0002 0.0001 -0.0001 0.0001 0.0000 0.0000 
R3 0.0004 -0.0011 0.0101 -0.0011 0.0004 -0.0002 0.0001 -0.0001 0.0001 0.0000 
R4 -0.0002 0.0004 -0.0011 0.0101 -0.0011 0.0004 -0.0002 0.0001 -0.0001 0.0001 
RS 0.0001 -0.0002 0.0004 -0.0011 0.0101 -0.0011 0.0004 -0.0002 0.0001 -0.0001 
R6 -0.0001 0.0001 -0.0002 0.0004 -0.0011 0.0101 -0.0011 0.0004 -0.0002 0.0001 
R7 0.0001 -0.0001 0.0001 -0.0002 0.0004 -0.0011 0.0101 -0.0011 0.0004 -0.0002 
R8 0.0000 0.0001 -0.0001 0.0001 -0.0002 0.0004 -0.0011 0.0101 -0.0011 0.0004 
R9 0.0000 0.0000 0.0001 -0.0001 0.0001 -0.0002 0.0004 -0.0011 0.0101 -0.0011 
R10 0.0000 0.0000 0.0000 0.0001 0.0001 0.0001 0.0002 0.0004 -0.0011 0.0100 


N2 
Table 4. Computer generated few values of Ay = vi pz1v, and A, = Pie 
S FINS. 
Iteration Iteration Iteration [Iteration Iteration [Iteration Iteration Iteration Iteration 
93 94 95 96 97 98 99 100 101 
Ay = vil p ws 78.05 66.12 56.37 48.73 42.98 38.86 36.10 34.53 34.02 
na N? 42.45 34.73 29.48 25.89 23.45 21.82 20.77 20.19 20.00 
s Vs Pnvs 
2.89 2.92 2.95 2.98 3.02 3.05 3.08 3.11 3.14 


Figure 9(a) illustrates high surface noise array gain (dB) vs @, for several S,(w,)/aZ levels of a 
MVDR Beamformer. Considering the same model except Sp, (w) interference comprises a white component 
with a spectral matrix o2Z] and a high-layer noise environment with a spatial spectral matrix mathematically 
represented as in (21): 


Sn (@) = So(w) [sincCkoap) 


+a (ax: — 1) sinc(koAp) 


- scos(koāp)} petan] 


(21) 

To get the optimal array gain and the traditional array gain, we repeat the same steps. Figure 9(b) illustrates 
high-level noise environment array gain (dB) vs 0, for several S,(w,)/o2Z levels of a MVDR Beamformer 
[39], [40]. 

In terms of hardware resources, speed (weight updating rate), algorithms implemented and 
applications proposed, the present work is compared with earlier recorded works [41], [42]. Table 5 shows 
the comparison. Although the proposed MVDR beamformer has high array gain for high surface noise 
compared to conventional. Comparing the results of the proposed study is not a reasonable estimate because 
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the literature research work was performed on multiple technologies (FPGA) and platforms, but the authors 
work are based on higher surface noise array gains than in the previous work [43]. The proposed work 
mainly focuses on improving the array gain throughput. According to the author’s knowledge, this is the first 
reported processor based MVDR architecture that enables real-time implementation of MVDR beamformer 
to achieve high array gain with high surface noise. Compared to the latest technology, it is novel in that it has 
the best array gain for high surface noise and high-level noise throughput. 


2 7 11.5 r J 
— MVDR 
= = Conventional 


Array gain (dB) 
Array gain (dB) 


— MVDR 


= = Conventional 


0 20 40 60 80 100 120 140 160 180 0 20 40 60 80 100 120 140 160 180 
4, (degrees) 6, (degrees) 


(a) (b) 


Figure 9. Array gain (dB) vs 8, for several S,,(w,)/a% levels of a MVDR: (a) beamformer for high surface 
noise and (b) beamformer for high layer noise 


Table 5. Performance comparison of MVDR with previous literature works 


Scheme No of Adaptive weight Updating Rate Clock Array Algorithm 
channels/Matrix computation (AWC) (K/Sec) cycles Gain in 
size /Matrix inversion (MI) dB 

Proposed 10x10 MI 120.5 1600 21.7 Processor based 
Kidav et al. [44] 32 AWC 132.2 1600 - Modified complex 
Ma et al. [45] 4x4 MI 17.62 - - Modified squared 

MInv/Sec givens rotation 

Rosado et al. [46] 23x23 MI 3.9 3441 - Qrd, modified gram- 

schmidt method 

Wang et al. [47] 23x23 MI - - 20.0 Maximum Diffuse 


noise gain beamformer 


5. CONCLUSION 

To this end, we propose an efficient optimal beamformer for discrete interference. It is robust not 
only to the volatility of the covariance matrix but also to vector steering. Based on the spatial spectrum 
distribution of the quasi-signal free environment, the interference plus noise covariance matrix can be 
reconstructed with the estimated DOA of the signal of interest. Based on the reconstructed covariance matrix, 
the estimated steering vector of the signal can be modified to optimize the output power of the array. Unlike 
other algorithms, the proposed optimal beamformer for discrete interference does not require standard 
steering vector constraints, thus greatly expanding its coverage. Simulation results show that the performance 
of the proposed optimal beamformer for discrete interference provides a significant performance 
improvement for large INRs in the sidelobe region and provides an acceptable beam pattern. MVDR 
beamformer produces solutions that distort the main lobe of the HPBW region and are prone to mismatch 
modeling. Considering multiple perturbation effects with numerous factors and enhancing MVDR array 
benefit for complex scenario processing, which is our future research work. 
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